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The double-resonant laser excitation of nitric oxide, cooled to 1 K in a seeded supersonic molecular 
beam, yields a gas of ^10 12 molecules cm -3 in a single selected Ryberg state. This population 
evolves to produce prompt free electrons and a durable cold plasma of electrons and intact NO + 
ions. This plasma travels with the molecular beam through a field free region to encounter a grid. 
The atomic weight of the expansion gas controls the beam velocity and hence the flight time from 
the interaction region to the grid. Monitoring electron production as the plasma traverses this grid 
measures its longitudinal width as a function of flight time. Comparing these widths to the width 
of the laser beam that defines the initial size of the illuminated volume allows us to gauge the rate 
of expansion of the plasma. We find that the plasma created from the evolution of a Rydberg gas 
of NO expands at a small but measurable rate, and that this rate of expansion accords with the 
Vlasov equations for an initial electron temperature of T e « 8 K. 



INTRODUCTION 

When a tunable pulsed laser promotes a substantial 
fraction of the atoms held in a magneto-optical trap 
(MOT) to a high Rydberg state (n > 30), the ensemble 
spontaneously evolves on a microsecond timescale to pro- 
duce an ultracold plasma [TJ [2j [3j HJ IS] • This phenomenon 
and the plasma that it forms have attracted a great deal 
of recent attention [6j [3 |8] , because these rarified labo- 
ratory systems can display a degree of charged-particle 
correlation otherwise found only in thermonuclear explo- 
sions and the cores of dense stars. 

Coulomb interactions act to govern the dynamics of 
the charged particles in a plasma when their electrostatic 
repulsion exceeds their thermal translational energy. The 
point at which this occurs depends on the density and 
temperature of the plasma [9]. The parameter, T, gauges 
the degree of correlation in terms of the dimensionless 
ratio, 

q 2 /47re a 

kT ' [ ' 

where q is the charge, and a, the Wigner-Seitz radius, 
relates to the particle density, p, by, 

4/3 Tra 3 = 1/p. (2) 

Conditions under which T substantially exceeds 1 can 
give rise to liquid-like or solid-like spatial correlations, 
leading ultimately to Coulomb crystallization and a 
quantum-state-detailed variant of the Mott insulator-to- 
metal phase transition [TO] . 

Ultracold plasmas formed in MOTs at densities of the 
order of 10 10 atoms cm -3 have an initial Wigner-Seitz 
radius of a = 3 jim. This average inter-particle spacing 
grows with a plasma expansion at rates that accord with 
electron temperatures falling from tens of K to the order 



of 10 K [TT]. At its maximum sometime during the 40 jis 
observation window, these conditions give rise to electron 
correlation parameters in the range of Y e « 0.1. 

Recently, we have demonstrated the evolution to an 
ultracold plasma by NO molecules entrained in a super- 
sonic expansion [12]. This development holds interest- 
ing potential for the study of correlation in mesoscopic 
charged-particle ensembles. The use of supersonic expan- 
sions substantially broadens the scope of systems that 
can be studied; any substance that can be volatilized — 
including large molecules and refractory metal clusters — 
can be entrained in a seeded supersonic molecular beam. 
Heavy-particle internal degrees of freedom present new 
dimensions for plasma relaxation dynamics, including an 
accommodation of electron energy in ion degrees of free- 
dom extending beyond the translational heat capacity 
to include rotational, and, in some cases vibrational and 
vibronic states. 

In our experiment, a seeded supersonic expansion pro- 
vides NO at a local density of 5 x 10 13 molecules cm -3 
and a moving-frame translational temperature less than 
1 K. Double-resonant laser excitation promotes about 
10 percent of these molecules to a selected high-Rydberg 
state, from which they evolve to form a remarkably 
durable ultracold plasma. At the Wigner-Seitz radius 
corresponding to this density (a = 500 nm), an electron 
temperature of 10 K would yield a correlation, T e w 4. 

Our initial observations provide no direct means to 
quantify the ion or electron temperatures in the fully- 
formed plasma. We do know, however, that the plasma 
expands very little during the 9 jis observational window 
of those experiments. 

We have now more than tripled this observation time. 
Replacing helium with neon, argon and krypton as ex- 
pansion gases into which we seed NO, we have reduced 
the beam velocity to substantially increase the flight time 
in the field- free region of our apparatus. By monitoring 
electron production as the plasma traverses the grid that 
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FIG. 1: Schematic diagram showing the path of the molecular 
beam from a differentially pumped source chamber to the 
interaction region through a skimmer to enter a system of 
three grids ending in a flight tube capped by a microchannel 
plate (MCP) detector. Drawn to approximate scale with a 
d to G2 spacing of 2 cm. 



defines the end of the field-free region, we measure the 
longitudinal width of this distribution of NO + ions and 
electrons as a function of flight times from 9 to 30 /as. 

Comparing these widths to the width of the laser beam 
that defines the initial size of the illuminated volume, we 
are now much better able to gauge the rate of expansion 
of the plasma. We find that the plasma created in a 
molecular beam from a Rydberg gas of NO does expand 
at a small but measurable rate, and that this rate of 
expansion accords with the Vlasov equations for an initial 
electron temperature of T e « 8 K. 



EXPERIMENTAL 

A pulsed jet of NO, seeded nominally at 10 percent in 
He, Ne, Ar or Kr at a backing pressure of 5 atm, expands 
through the 0.5 mm diameter nozzle of a Series 9 pulse 
valve (Parker Hannifin). Two centimeters downstream, 
the free jet passes through an electroformed Ni skimmer 
to enter a mu-metal shielded electron spectrometer as a 1 
mm diameter differentially pumped supersonic molecular 
beam. 

Figure [I] diagrams our two-stage electron spectrom- 
eter. In the first stage, a pair of grounded grids, Gi 
and G2, held perpendicular to the axis of the molec- 
ular beam, define a field-free laser-molecular-beam in- 
teraction region. Tunable pulses from two, frequency- 
doubled Nd:YAG-pumped dye lasers overlap to intersect 
the molecular beam halfway between the plates. 

Here, NO molecules absorb light in two resonant steps 
to reach a selected high Rydberg state. These high- 
Rydberg NO molecules interact to form a plasma. This 
plasma volume moves with the velocity of the molecu- 
lar beam to traverse G2, at which point it encounters 
a 60 V cm -1 field. This field gradient extracts and ac- 
celerates plasma electrons to a multichannel plate de- 
tector. On the timescale of our measurement, electrons 
extracted at G2 appear instantaneously at the detector, 
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FIG. 2: Excitation scheme for the population of nf Rydberg 
states by double-resonant transitions via the 2 £ + ,v = 0,iV = 



0, J = |(+) state of NO. 



and the profile of this signal thus integrates transverse 
slices of the electron density distribution of the plasma 
as it passes through G2. 

Figure |2] shows an energy diagram for the excitation of 
NO. An initial laser pulse (coi) pumps the transition from 
the ground, X 2 n state to rotational substates of the the 
v = vibrational level of the first electronically excited 
state, 2 S + . A second laser (^2), timed to coincide with 
uj\ excites from a chosen level in this gateway system to 
a selected high Rydberg state situated below the lowest 
ionization threshold. 

At high laser pulse energies two-photon absorption of 
U 1 alone ionizes NO. Monitoring the electron signal as 
a function of the uj\ wavelength under these conditions 
produces a rotationally resolved resonant ionization spec- 
trum of the v = band of the 2 S + state. We analyze 
the intensities of lines in this spectrum to estimate the 
rotational temperature of NO under our expansion con- 
ditions. 

For plasma experiments, we attenuate uj\ using a pair 
of Glan- Taylor polarizing prisms to eliminate any back- 
ground signal from uo\ alone, and then counter propagate 
unfocussed uo\ and CJ2 laser beams to intersect the molec- 
ular beam as indicated in Figure [I] This illumination 
geometry crosses a cylindrical photon field with a col- 
limated molecular beam to produce a prolate ellipsoid 
excitation volume that propagates sideways with the ve- 
locity of the beam toward G2. The diameter of uo\ at 
its intersection with the molecular beam determines the 
initial equatorial width of this ellipsoid. 

We experimentally gauge this limiting laser spot size 
by razor-blade tomography. Figure [3] traces the inte- 
grated intensity of uj\ as a function of the position of 
a razor blade scanned through the beam waist at a dis- 
tance corresponding to the point at which the laser pulse 
crosses the axis of the molecular beam. This measure- 
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FIG. 3: Horizontal profile of u)\ the smaller of the two laser 
beam diameters, which determines the width of the plasma 
in the direction of molecular beam propagation . 



merit yields a value of 758 /im for the initial width of 
the active excitation volume along the propagation axis 
of the beam. The divergence of the molecular beam af- 
ter the skimmer determines a cross-beam length for this 
illuminated cylinder of approximately 5 mm. 
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FIG. 4: Resonant 1+1 ionization spectrum of the X 2 II v" — 
0) to A 2 S + v — 0) transition in NO. This pattern of intensi- 
ties fits a ground state population distribution corresponding 
to a rotational temperature, Tr = 2.5 K. 



RESULTS 

A seeded supersonic expansion under the present con- 
ditions of backing pressure nozzle diameter cools a di- 
atomic gas to a parallel temperature, Tjj 50 , of about 
700 mK, where Tjj 50 represents the second moment of 
the velocity distribution function along the axis of the 
molecular beam [13] . Supersonic expansion cools molec- 
ular rotation as well, but the cross section for rotational 
relaxation is smaller than that for translation. As a re- 
sult, when collisions cease, the terminal rotational tem- 
perature, TgP, exceeds . Figure |4j shows a uo\ scan of 
the resonant ionization spectrum of the X 2 II v" = 0) to 
A 2 S + v' — 0) transition of NO under the conditions of 
our experiment, from which we estimate = 2.bK . 

From the nozzle diameter and backing pressure, we es- 
timate the centerline density of particles in the molecular 
beam to be 4.8 x 10 14 cm" 3 QUE . Tne seedm g ratio of 
NO is 0.1. At a rotational temperature of 3 K, 87 percent 
of these molecules populate the two parity components 
of the rotational ground state. Two saturated steps of 
laser excitation transfer 12.5 percent of this population 
to a parity selected high Rydberg state. In combination, 
these factors predict the density of excited NO available 
for plasma formation to be 5 x 10 12 cm -3 . 

To prepare a plasma, we reduce the power of uo\ , set its 
frequency on the Qi(^) line of this system (X 2 II v" = 
0, N" = 1, J" = §(-) to 2 S+ v' = 0,N f = 0, J' = §(+)), 
and tune UO2 so that the total uj\ + 002 energy approaches 
the adiabatic ionization threshold of NO. Figure [5] shows 
electron signal waveforms obtained for expansions of NO 
in He, Ne and Ar, consisting of a prompt electron signal 
that appears within a few nanoseconds of 002 followed 
tens of microseconds later by a broader late signal. 

In previous work, we have established that these wave- 
forms reflect the production of an ultracold plasma of 



electrons and molecular nitric oxide cations [12] . The 
prompt signal arises from electrons released early in the 
formation of the plasma. The late signal marks the pas- 
sage of the quasi-neutral plasma through the grid, G2 at 
the laboratory velocity of the molecular beam. 

Gating on the late signal and scanning the wavelength 
of UO2 1 we obtain a spectrum of closely spaced resonances, 
which is shown for a backing gas of Ar in Figure [6] We 
assign this spectrum to lines in the nf{2) Rydberg series 
[T5] . where the number 2 refers to the rotational level 
of NO + to which the series converges. The time traces 
shown in Figure [5] are obtained for 002 tuned to 52/(2). 

For each expansion gas, the measured arrival time of 
the plasma at G2, coupled with the known flight distance 
from the laser interaction region, yields a sharply defined 
laboratory-frame average velocitiy, < v z > by which to 
convert the observed temporal width to a spatial width 
along the propagation axis of the molecular beam. Fig- 
ure [7] plots the late-peak waveforms from Figure [5] trans- 
formed to distance. We interpret these growing widths to 
represent the plasma density distribution as it expands 
in the spatial dimensions transverse to the laser illumi- 
nation axis. Distances measured in the z-dimension fit 
well with a representation of this expansion in terms of 
a two-dimensional Gaussian: 



n ^ z ^ t) = ^m exp[ -^w~ ] (3) 

for which V — 2\^2ln2a(t) defines the full- width at half- 
maximum of the plasma density distribution in space, 
as defined by the time- varying standard deviation, cr(t). 
Table [T] lists the arrival times and spatial widths of the 
electron-density images obtained in this way for all four 
expansion gases used in this study. 
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FIG. 5: Electron signal waveforms appearing at G2 follow- 
ing the substantial promotion of NO molecules to the 52/(2) 
Rydberg state in seeded supersonic expansions of He, Ne, Ar 
and Kr (from bottom to top) using the experimental flight 
path diagrammed in Figure^ 
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FIG. 6: Resonances in the late signal observed scanning cj 2 - 
Spectrum assigned to lines in the nf Rydberg series con- 
verging to the rotational level, N + = 2 in the vibrational 
ground state of NO + X 52/(2) falls at an 002 energy of 

30,493.9 cm -1 , with a convergence limit of 30,534.5 cm -1 . 



DISCUSSION 

By choice of the expansion gas, we systematically vary 
the laboratory velocity of nitric oxide in a differentially 
pumped seeded supersonic molecular beam. Crossing 
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FIG. 7: Late electron-signal temporal waveforms from Fig- 
ure [5] transformed to yield spatial plasma distributions, dis- 
played according to peak arrival time at G2. Solid curves 
show Gaussian fits with residuals below. The bottom figure 
is derived from the measurement of laser width shown in Fig- 
ure [3] above which are plotted data for He, Ne. Ar and Kr 
carrier gases. 
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TABLE I: Arrival times, estimated velocities and widths ob- 
tained by fitting Gaussian distributions to the spatial wave- 
forms pictured in Figure [7] as transformed from the flight- 
time dependent electron signals shown in Figure [5] 



this beam with the output of two dye lasers, we use 
double-resonant excitation to transform a prolate volume 
element with a minor axis diameter smaller than 1 mm 
into a dense distribution of NO molecules excited to a se- 
lected high-Rydberg state substantially below the lowest 
ionization threshold. 

This population of high-Rydberg molecules evolves to 
form a plasma of NO + cations and electrons. The volume 
element containing this plasma moves through a field- 
free region with the propagation velocity of the molec- 
ular beam to traverse a perpendicular grid and enter a 
60 V/cm field gradient. Transmission of this volume ele- 
ment through the grid disrupts the plasma to produce a 
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FIG. 8: Plot of plasma FWHM at G 2 as a function of flight 
time, compared with width predicted after 10 /us for ambipolar 
expansion at T e = 40 K (open dot) and Vlasov expansion of 
a spherical Gaussian plasma for T e (0) = 40 K (bold line). 
The line through the data represents a Vlasov fit for T e (0) 
= 7.8 K. Error bars reflect the experimental uncertainty in 
the measured flight distance to G2, which exceeds the error 
introduced by using velocity averages to transform time to 
distance. 



waveform of promptly detected electrons that maps the 
plasma charge-density distribution along the axis of the 
molecular beam. 

Converting these temporal images to displacements in 
space, we measure the physical width of the plasma along 
its laboratory propagation axis as a function of flight time 
from 9 to 30 /is. Figure [8] plots these results. With ref- 
erence to the diameter of the initial illumination volume, 
we find that in 30 /is this measure of the plasma width 
approximately doubles. 

The expansion rate of an ultracold plasma relates prin- 
cipally to the thermal energy of its electrons. Consider- 
able effort in the field of atomic ultracold plasma research 
has focused on the determination of the electron temper- 
ature over the life of the plasma and its relation to plasma 
expansion [TTJ [16] [T7] [18] [19] . These measurements have 
generally supported the idea that, for plasmas prepared 
above threshold, the quantity, T e = 2E e /3kB, reason- 
ably describes the initial electron temperature, so long 
as E e , the excess energy, exceeds ~ 50 cm -1 . For plas- 
mas produced by photoexcitation tuned lower in energy 
to release free electrons with less translational energy, 
and lower still to populate high-Rydberg states, electron 
heating mechanisms dominate to maintain initial electron 
temperatures in the range of 40 K. 

Electrons released by plasma formation expand much 
more rapidly than the colder (T^ « 1 K), heavier pos- 
itive ions. However, the attractive force of a develop- 
ing space charge acts to constrain the expanding electron 
spatial distribution to grow no faster than the timescale 
for cation motion. By this Coulomb coupling, the ther- 
mal energy of the electrons drives overall plasma expan- 
sion, viz: 



^M iV 2 = h B T e 



(4) 



For NO + in a plasma with an electron temperature of 
40 K, this expression predicts an expansion rate of 100 
/im /is -1 . An open dot on Figure [8] shows the width after 
10 /is predicted by Eq. [4] for a plasma with an initial 
diameter, <r(0) = 758 /im. Clearly, the size expected 
for ambipolar expansion at T e (0) = 40 K substantially 
exceeds our observations. But, we should expect this 
simple picture to overestimate the expansion rate. By 
expanding, the electron charge distribution does work on 
the cations, so the temperature of the system, and thus 
the force driving its expansion, should fall with time. 

A fuller account of the coupled evolution of the parti- 
cle density and energy distribution functions is provided 
by the Vlasov equations [20 , which form the founda- 
tion of the kinetic theory of plasmas, and have analytical 
forms for many types of low-density collisionless plasmas 
[2T] . Laha and coworkers [11 have applied an analytical 
self-similar solution to describe the expansion of a quasi- 
neutral ultracold plasma configured to have a spherically 
symmetric, Gaussian ion-density distribution. The evo- 
lution timescale for such a plasma is given by: 



m iC r(0) 2 



MT e (0)+Ti(0)] ' 



(5) 



where the quantities denoted T a (0) refer to the initial 
temperature of the electrons, T e (0), and the ions, 1^(0). 
mi is the mass of the ions. 

The electron and ion temperatures fall with time ac- 
cording to 



T a (t) 



T a (0) 



(6) 



and the diameter of the plasma expands as 



a(t) = Ja(0y 



72~ 



(7) 



This simple formalism conforms very well with experi- 
mental measurements of ion positions and correlated ex- 
pansion velocities in quasi-neutral ultracold plasmas of 
strontium prepared to have electron temperatures in the 
range of T e (0) = 40 K and above [EE]. 

For reference, the bold line on Figure [8] gives the time- 
dependent width predicted by Eq.[7]for a plasma with the 
ion mass of NO, an initial width of 785 /im and initial 
electron and ion temperatures of T e (0) = 40 K and T^(0) 
= 1 K. The curve shown significantly overestimates the 
rate at which our plasma expands. 
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The collimated laser-crossed molecular-beam excita- 
tion geometry of our experiment produces a prolate ellip- 
soidal plasma volume, as opposed to a Gaussian sphere. 
Nevertheless, tomography shows that our laser inten- 
sity distribution in the yz plane is Gaussian, and slice 
measurements of the electron signal along the molecular 
beam propagation axis show that this distribution of ex- 
citation intensity gives rise to a Gaussian plasma density 
distribution in the equatorial dimension. 

Experiments by Cummings et al. [5 , form an ellipsoidal 
plasma volume of similar aspect ratio in a calcium MOT. 
They determine an equatorial expansion rate by model- 
ing the time-dependent Ca + laser-induced fluorescence 
signal. Though the analytical self-similar solution of the 
Vlasov equations for a spherically symmetric Gaussian 
plasma do not extend to this geometry, the expansion 
rate deduced by Cummings et al. [5 for calcium accords 
almost exactly with that found experimentally for the 
spherical strontium plasma of Laha et al. [11 when 
matched for initial electron energy and scaled to account 
for the ion mass difference. 

This suggests that a simple fit of the spherical formal- 
ism to the rate of expansion in the equatorial dimension 
provides a reasonable gauge of electron temperature for 
Gaussian ellipsoidal plasma geometries. 

The line drawn through the data on figure [8] represents 
a Vlasov fit to the equatorial plasma diameters measured 
in our experiment, returning initial electron temperature 
of T e (0) = 7.8 K (for an assumed ion temperature of 
7i(0) = 1 K). Remarkably, by this fit, the Vlasov equa- 
tions hold that after 30 /is the electron temperature of our 
plasma has fallen to 1.6 K. By this point, plasma expan- 
sion has doubled the Wigner-Seitz radius, but the pre- 
dicted electron correlation remains substantial, Y e ^10. 

This plasma forms following nitric oxide Rydberg- 
Rydberg Penning ionization collisions. Some initially 
formed electrons leave the excitation volume immedi- 
ately, creating an electrostatic trap in which an electron- 
NO* avalanche ionizes the remaining Rydberg molecules. 
Pillet, Pohl and coworkers have recognized the collisional 
ionization of Rydberg atoms as a cooling mechanism for 
electrons in ultracold plasmas [22]. MOT experiments 
appear to show that the properties of atomic ultracold 
plasmas vary continuously with excitation energy tuned 
across the ionization threshold (c.f. Figure 5 in reference 
[5]). Thus, we expect the plasma formed from Rydberg 
NO molecules to exhibit a relatively low electron temper- 
ature. 

However, as mentioned above, efficient heating mecha- 
nisms appear to elevate electron temperatures, even in 
systems prepared very near threshold, to levels much 
higher than seem to be evident here. Chief among these 
for the lowest initial values of T e is three-body recombi- 
nation [T8] . 

In simple physical terms, three-body recombination oc- 
curs when an electron approaches an ion to within a crit- 



ical radius, ry, during which time a collision with a sec- 
ond electron carries away the binding energy of the pair. 
By this mechanism, any ion-electron collision resonance 
bound at least by the energy to which the electron colli- 
sionally thermalizes can remain bound. Thus, T e deter- 
mines the largest orbital radius stabilized by three-body 
recombination (the Thomson radius [23]) by the energy 
balance: 



\k B T e 



(8) 



This largest, lowest binding energy orbital dominates 
recombination by virtue of the cross sections for elec- 
tron capture and collisional thermalization, which scale 
as r 2 . Thus, the rate of three-body recombination vastly 
accelerates as the electron temperature approaches K. 
Various collision theory formulations developed to model 
three-body recombination consistently yield third order 
rate constants that scale with temperature as T~ 4,5 
US [23] El EU EH [27] . Calculations for ultracold atomic 
plasmas show that three-body recombination makes it 
difficult to form a plasma with an electron temperature 
less than 25 K [T9], 

The very high rates predicted for low-temperature 
three-body recombination follow directly from the very 
large values of Tt sampled in theory by low-energy re- 
combining and deactivating electrons, or equivalently, the 
very high-values of maximum principal quantum num- 
ber, n max (and accompanying orbital degeneracy), set 
by temperature in detailed-balance integrals. 

However, conventional formulations of the three-body 
recombination rate neglect the effect of plasma density on 
the accessible scale of tt or n max . For an electron tem- 
perature, T e = 2 K, the average kinetic energy, 2.6 meV, 
dictates a Thomson radius of 560 nm. For a typical MOT 



plasma with a density of 3 x 10 



15 



this lies well 



within the Wigner-Seitz radius of 4 /im. By contrast, we 
estimate the plasma formed under our molecular beam 
conditions to have an ion density of 2 x 10 18 m -3 , or a 
Wigner-Seitz radius of 500 nm. Models incorporating the 
effects of ion density on high-n capture yield predicted 
rate coefficients that scale as T e -1 , significantly reducing 
the predicted rate of three-body recombination [28] [29] . 

Thus, it appears that conditions are right in our beam 
experiment to suppress three-body recombination and its 
attendant electron heating effects, diminishing the rate 
of plasma expansion in comparison with atomic systems 
prepared in MOTs. 

A suppression of ion-electron recombination seems also 
apparent in the absence of dissociative recombination. 
We have noted previously [12 that a plasma lifetime of 
9 [is exceeds expectations based on known cross sections 
for the dissociative recombination of NO + [30 . Present 
results include NO + plasmas that maintain a high degree 
of integrity for as long as 30 /is. 
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